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Abstract 

We have studied the entropy- driven mechanism leading to stationary patterns for- 
mation in stochastic systems with local dynamics and non-Fickian diffusion. We 
have shown that a multiplicative noise fulfilling a fluctuation-dissipation relation 
is able to induce and sustain stationary structures with its intensity growth. It 
was found that at small and large noise intensities the system is characterized by 
unstable homogeneous states. At intermediate values of the noise intensity three 
types of patterns are possible: nucleation, spinodal decomposition and stripes with 
liner defects (dislocations). Our analytical investigations are verified by computer 
simulations. 
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1 Introduction 

It is well known that the interaction of noise and nonlinearity attracts a lot 
of interest in nonequilibrium systems theory [1]. Most of articles, concerning 
above noise induced phenomena, relate to problems of external noise influence: 
noise induced transitions in zero-dimensional systems [2] ; noise-induced phase 
transitions [3,4]; stochastic resonance [5]; noise sustained patterns [6-9], etc. 
Analytically, numerically and in natural experiments it was found that an ex- 
ternal noise plays an organizing role only if its amplitude depends on the field 
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variable — multiplicative noise (see, for example, Refs. [2,4,10,11]). If more 
than one noise source is introduced into the system dynamics, then the noise 
cross-correlations start to play a crucial role in ordering processes: a kind of 
phase transitions principally changes [12,13]. From a fundamental viewpoint 
such effects has a dynamical origin: in a short-time limit external fluctuations 
destabilize a disordered homogeneous state. An analytic description of ex- 
tended systems with external fluctuations is provided with an approximately 
known stationary distribution function. 

The possibility of noise induced patterns formation has been studied in last 
two decades (see, for example, Refs. [4,6-9,14-17]). One of the most interesting 
problem when the noise sustains patterns is widely discussed recently [17-21]. 
In the most of works, concerning this problem, external fluctuations influ- 
ence on the dynamics of systems is considered. Recently a new mechanism for 
noise induced patterns formation was discovered. Its origin is in a relaxational 
dynamics governed by the field-dependent mobility/kinetic coefficient [22]. It 
was shown that corresponding self-organization processes are not related to 
a short-time instability of a homogeneous/disordered state. Here the princi- 
ple role is played by entropy variations caused by the field- dependent mobil- 
ity [23-28]. Thus, an ordered behaviour appears due to the balance between 
relaxing forces moving the system to the homogeneous state, and fluctuations 
dependent on the field variable pulling the system away from the disordered 
state. This mechanism belongs to a kind of entropy driven phase transitions 
that is an extension of noise induced unimodal-bimodal transitions in zero- 
dimensional systems [2]. In such systems a dissipation is related to fluctua- 
tions according the fluctuation dissipation relation. Hence, one gets an internal 
multiplicative noise with intensity reduced to the bath temperature. For such 
class of stochastic systems the corresponding distribution function, free en- 
ergy and associated effective potential are known exactly. Therefore, coherent 
structures/patterns appeared in the course of self-organization processes can 
be analyzed without any dynamical reference. 

In this paper we will study a behaviour of extended stochastic systems when 
the noise can induce or sustain spatial patterns according to entropy driven 
mechanism. To this end we consider a model of reaction-diffusion system of the 
field-dependent mobility with non-Fickian diffusion and internal fluctuations 
of the multiplicative character. Reaction-diffusion systems play an important 
role in the study of generic spatiotemporal behaviour of far from equilibrium 
systems. These models are usually applied to describe an extended systems 
when the physical space is divided in to the small cells and dynamics of each 
cell is described by a variable type of concentration. The model to be consid- 
ered concerns the local dynamics of the variables inside each cell and the trans- 
port phenomena between cells. The local dynamics is determined by so-called 
chemical reactions inside each cell. The transport is caused by non-Fickian dif- 
fusion. We will study stationary spatial inhomogeneities with lengths bigger 
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than the linear dimensions of the cells in such system. 

Formally, reaction-diffusion models are described by the equation for the field 
x = x(r,t) in the form d t x = f(x) — VJ, where f(x) stands for the lo- 
cal dynamics, J is the flux for transport phenomena, V = d r . Considering 
the non-Fickian diffusion, we exploit a gradient of molecular interactions po- 
tential U(r) (U(r) = — J u(r — r / )x(r / )dr / , where the spherically symmetric 
interaction potential u(r) between molecules separated by a distance |r| is in- 
troduced). The corresponding force given by the gradient of U(r) governs the 
transport phenomena. In the case of the small interaction radius compared to 
the diffusion length the concentration field x will not vary significantly within 
the interaction radius. It allows to approximate the integral by kx + /3V 2 x, 
where k = f w(r)dr, (3 = (1/2) / |r| 2 w(r)dr (j3 ~ «r 2 /2, where r c is a correla- 
tion radius); Jrw(r)dr = due to symmetrical properties of the interaction 
potential [29]. Therefore, the obtained flux allows to describe phase separation 
processes with mutual (lateral) interactions; the combined model with local 
dynamics can be used to consider the spatial patterns induced by the fluctu- 
ations of the bath. With the help of such formalism one can find an explicit 
form of the stationary distribution functional and study spatial patterns using 
a variational principle. Analytical results we verify by computer simulations. 

The paper is organized in the following manner. In the Section II we intro- 
duce a general model for stochastic reaction-diffusion systems and present the 
theoretical approach is useful in our consideration. In Section III we discuss 
the model and assumptions related to our study. Section IV deals with a the- 
oretical investigation of stationary noise induced/sustained patterns. Finally 
we summarize our results in Section V. 



2 Main equations 

Let us start with the typical deterministic equation for the reaction-diffusion 
models in the form 



where x stands collectively for the variables, f[x] for the local dynamics inside 
each cell, J[x] is the flux proportional to the conjugate thermodynamic force, 
arising from the spatial variation of the chemical potential /j,[x(r,t)]. One can 
exploit a definition J[x(r,t)] = — D e f[x(r, t)] Vfi[x(r, t)], where the chemical 
potential fi = 5J-/5x is introduced through the free energy J-[x]; D e f[x] is an 
effective diffusion coefficient related to the field- dependent mobility. In such a 
case, the deterministic evolution equation for the field x reads 



d t x = f[x(r,t)]-VJ[x(r,t)}, 




d t x = f[x]+V(D ef [x]V f i[x]). 



(2) 
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Formally, Eq.(2) can be written in a variational form as 

° tX ~ D e ,\x] Sx • (6> 

where the potential functional U[x\ is represented through the definitions for 
f[x], D e f[x] and T[x\. It plays a role of a Lyapunov functional for the deter- 
ministic dynamics. An explicit form for the functional U[x] can be found only 
if fi = x. Indeed, here one has 

U[x] = J dr j- J f[x']D ef [x'}dx' + (D ef [x]Vxf/2 

Generally, if the chemical potential ji is a function of both x and Vx, then 
one can obtain only its first variation, 

SU[x] = - J dr5x {f[x]D ef [x} + D ef [x]V (D ef [x] V/i[x])} , (4) 

which after substitution into Eq.(3) gives Eq.(2) immediately. 

Considering the system under real conditions of thermal bath influence, one 
has to take into account corresponding fluctuations. In stochastic analysis we 
introduce a related multiplicative noise in an ad hoc form, following Ref. [17], 



where the fluctuation dissipation relation holds, i.e. (£(r, t)£(r', t')) = 2a 2 5(r — 
r')5(t — t')] (£(r, t)) = 0; a 2 is the intensity of a Gaussian noise ^. 

Using the Stratonovich interpretation of the stochastic equation (5), the sta- 
tionary solution of the corresponding Fokker-Planck equation can be written 
as follows [23]: 

Vst[x] oc exp f-^P) , (6) 
where the effective potential functional C/ e /[x] is of the form 

U ef [x] = U[x] - y J drlnZ? e/ [x]. (7) 

It follows, that the stationary distribution functional V st or the effective func- 
tional U e f are obtained exactly: the form of initial functional U is supposed 
to be known, the second term in U e f can be calculated if needed. Let us note, 
if we assume that U plays a role of an effective free energy functional, then 
rewriting the integral in Eq.(7) as S e f = J dr ln(L> e /[x]) _1 , the expression (7) 
can be transformed into the thermodynamic relation between free energy, in- 
ternal energy and entropy functionals: U e f = U + a 2 S e f. Therefore, according 
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to such relation the noise intensity a 2 addresses to an effective temperature of 
the bath, whereas S e f plays a role of the effective entropy. Such situation is 
well known in the stochastic systems theory. It appears when the multiplica- 
tive fluctuations corresponded to the internal noise. The last one results to 
the entropy change that yields entropy driven phase transitions [18,23,26,28]. 
In this paper we will not discuss above phase transitions, we will consider an 
ability of the noise to sustain or induce formation of spatial structures. 

Stationary noise sustained structures x s (r) correspond to the extreme of the 
functional U e f\x\. To find an associated equation for x s {y) we make the first 
variation of the effective functional with respect to x and equal it to zero, that 
is 



M = -/dr fa {| 



SU ef [x] = - I dr6x { [f[x]D ef [x] + D ef [x]V (D ef [x]V ^[x]) + £ 

(8) 

As a result, the corresponding equation takes the form 

V (D ef [x]V»[x]) + f[x] + 2 °' x] dDe d f } X] = 0. (9) 

ef L J 

Homogeneous states are defined as solutions of the reduced equation 

„ o a 2 dD P f , . 

fI * = -2-te- (10) 

From the mathematical viewpoint solutions of Eq.(10) give extreme positions 
of the function U e f(x) obtained under supposition that x(r) = const. 




3 The model 



In order to investigate a possibility of patterns formation we use the well 
known assumptions for the diffusion coefficient D e f. the diffusion coefficient 
decreases from the constant value with the field x growth. Following Ref. [22], 
the generalized formula for D e f can be written in the form D e f(x) = (1 — x 2 ) a , 
where a > 0. Under such an assumption the diffusion D e f decreases slightly 
in the vicinity of x = and goes to zero with an increase in x by absolute 
value. At small a one can use an approximation 

D ef =— — 2 , a>0. (11) 
1 + ax 2 

Such a construction assumes that fluctuations can not disappear with an in- 
crease in x, moreover some analytic solutions of the problem can be obtained. 
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The reaction term f(x) can be defined according to a chemical kinetics and 
generally is represented through a potential function, V{x) in the standard 
way: f(x) = —dV/dx. In this work we study a case of nonlinear potential 
force of the form f(x) = — Yli(x — x ^), where the set {x ^} corresponds to 
zero values of the force and relates to stationary points of the deterministic 
system. In our approach this force associates with the potential 



V(x) = y + ^ 3 " \*\ (12) 



here fi and e are constants that control the system dynamics. The potential 
(12) has three extrema placed at Xq = — fi/2 ± \J \i 2 + 4e/2, and at xq = 0. 
A spinodal is given by equation e = —fi 2 /4. A prototype model of chemical 
reactions is A ^ B with transient reactions A + 2X 3X , X ^ B, where 
first reaction occur with constant rates k\, ki, the second one is realized with 
ks, k^ [30]. Here x G [—1,1] measures deviations of spices X concentration 
from the constant fixed value controlled by parameters e and fi, which relate 
to the constant rates k\, k 2 , k 3 , k 4 . 

To make the quantitative analysis let us use the chemical potential in the 
standard form \i = bT '/5x, where the free energy functional is assumed in the 
Ginzburg-Landau form 



0(x) + \ (Vx) 2 



(13) 



4>{x)) the local potential; (3 > is an inhomogeneity constant. In oder to 
make a general description of patterns formation scenario let us consider the 
prototype model of phase separation where the local potential is of the form 



<t> = -- 2 x' + \x\ (14) 



where k is a control parameter. From the formal view point one can consider 
two different models. Indeed, in the simplest case of lateral interactions one 
can use the local potential in the form <f){x) = —fx 2 , k > [29]. A case of 
ordinary phase separation scenario corresponds to the x 4 -construction, given 
by Eq.(14). Next, we derive a formalism for a general model and compare 
results for these two models of phase separation. 
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4 Results 



4.1 Noise induced transitions 



Let us note that the noise influence leads to a short-time instability in our 
model. Indeed, the linear stability analysis of the linearized Langevin equation 
(5) allows to find an evolution equation for the structure function S(k., t) = 
(xk(t)x_k(t)}, where x(k,t) = J x(r, t)e lkr dr is the Fourier transform of the 
concentration field, (...) represents an ensemble average over noise. Consid- 
ering the state x = 0, a spherically averaged structure function S(k,t) = 
J S(k, t)dftk (f^ is the hyperspherical shell of radius k) evolves according to 
a linear equation 

1 t\R(h f\ 

= -co(k)S(k,t) + a 2 , (15) 



2 dt 

where cu(k) = k 2 ((3k 2 — k) — e — aa 2 . It follows that internal multiplicative 
noise leads to instability of the null state. The stationary value of the structure 
function is 



a 2 



Sst{k) ~ k 2 ((3k*-K)-s-aa 2 - (16) 
It is seen that at short-time scales with an increase in the noise intensity the 
peak position of the function S st (k) is shifted toward large values of the wave 
vector. 

To make an appropriate analysis of patterns formation scenario, let us con- 
sider a case of a homogeneous (zero-dimensional) system, assuming stochastic 
variable depending on the time only, i.e. x(r, t) = x(t). The main attention will 
be paid to stationary states investigation. In the case under consideration the 
stationary probability density function is of the form P st oc exp( — U e f(x) /a 2 ), 
where the functional is reduced to the function, U e f[x] — »■ U e f(x). Hence, the 
effective potential is U e f(x) = — J dx' f(x)D e f(x) — (cr 2 /2) In D e f(x). Inserting 
definitions for / and D e f into this construction, we obtain 

U e f(x) = —x 2 +—x — — (e+ — ln(l+ax 2 ) ^arctan(a\/x)+— ln(l+ax 2 ). 

2a a 2a V otJ a 6 ' 1 2 

(17) 

To consider a case x G [—1,1] let us assume values of both e and fi to locate 
a minimum U e f(x-) at x < 0, a minimum U e f(x+) we locate at x > 0. An 
appropriate choice is e = 0.2, \i = —0.5, a > 0; the corresponding dependence 
U e f(x,a 2 ) is shown in Fig.l. To find most probable states we solve the prob- 
lem dll e f(x)/dx = 0. The corresponding equation takes the form (10). As it 
follows from our consideration the functional dependence of D e f(x) leads to 
bifurcations (a number of extrema of the effective potential changing). Hence, 
solutions of Eq.(10) allow to find bifurcation diagram illustrating extrema po- 
sitions of the effective potential. In order to calculate critical points related 
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Fig. 1. The effective potential U e f(x,a 2 ) at a = 0.2, e = 0.2, \i = —0.5. 

to the bifurcations one needs to differentiate Eq.(10) with respect to x and 
after solve both obtained equation and Eq.(10) simultaneously. Solutions of 
Eq.(10) show that a root xq = exists always. Another two roots 

[L t 1 

are realized if a 2 < a 2 where 



x ± = -^±-^ 2 + Ae-Aa 2 a (18) 



-I 




(19) 



At a 2 = a 2 solutions x_ and x + degenerate, and at a 2 > a 2 only trivial one, 
xq = 0, remains. The corresponding dependencies x±(a 2 ) are shown in Fig. 2a. 
To understand transformations of the system states let us use the noise in- 
duced transitions formalism [2]. As it follows from naive consideration, the 
bimodal stationary distribution P st {x) oc exp(— U e f(x)/a 2 ) becomes unimodal 
with an increase in the noise intensity a 2 . In the case under consideration 
above transition occurs in the following manner. At a 2 = a form of the 
effective potential U e f is identical topologically to a form of the initial poten- 
tial V(x). With an increase in the noise intensity a 2 a minimum of U e f(x-) 
located at x_ tends to zero, at a 2 = a 2 = eja the effective potential has 
a double degenerated point, Xq — X _ = 0. Therefore, the values a 2 define a 
spinodal curve. At a 2 < a 2 < ctq the point xq relates to a minimum, whereas 
X- defines a maximum position of the function U e f. These two minima dif- 
fer in depth, i.e. U e f(0) > U e f(x+). At a 2 = ctq one has U e f(0) = U e f(x + ), 
therefore, o"q defines a coexistence line (binodal). With a further increase in 
a 2 we get U e f(0) < U e f(x + ). An equality U e f(xJ) = U e f(x + ) is satisfied at 
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Fig. 2. Bifurcation diagram (a) and phase diagram (b) for noise induced transitions 
(a change number of extrema of the function U e Ax)) at a = 0.2, e = 0.2, (jl = —0.5. 
In plot (a): solid lines define stable states, dashed line corresponds to unstable 
solution. In plot (b): solid line corresponds to a 2 values, dashed and dotted lines 
relate to a 2 and ctq, respectively. 

a 2 = a 2 , hence the bifurcation point, a 2 , define another spinodal. At a 2 > a 2 
the effective potential has one well only. Therefore, in such a noise induced 
transition we have shift of the potential extreme, transformation of the global 
minimum into a local one, loss of its stability and, finally, change a number 
of extreme of the function U e f. The corresponding diagram illustrating above 
situation is shown in Fig. 2b. 



4.2 Noise induced patterns 



Next, let us consider a possibility of the system to manifest stationary struc- 
tures formation. To this end we assume the spatial dependence of the stochas- 
tic field and solve the variation problem, 5U e f[x\/5x = 0. The corresponding 
equation takes the form of Eq.(9). Let us rewrite Eq.(9) in a more convenient 
form 



^-gv^-g(v^ 



dD ef d 2 (f) 2 a 2 8D ef 

dx dx 2{ } 1 2D 2 ef dx ' 



(20) 

As it follows from our consideration, the system stationary behaviour can be 
described in the 4-dimensional space (x, y, z, u) where y = Vx, z = V 



V 3 x. Using results form the homogeneous system analysis one can state that 
there are three fixed points describing homogeneous states with coordinates 
0,0,0), (0,0,0,0) and (x + , 0,0,0) at a 2 < a 2 and the unique fixed point 
(0,0,0,0) at a 2 > a 2 c . 

Firstly, let us investigate the stability of stationary solutions of Eq.(20) in 
the r-space. To this end one assumes x(r) to be in the form x(r) = x(0) + 
<5exp(Ar), 5 <C 1 is a small perturbation in the vicinity of the corresponding 
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homogeneous solution. Inserting this assumption with corresponding expres- 
sions for all possible derivatives into Eq.(20), one gets the equation 



= (21) 

for eigenvalues A = 7±ifc of the Jacobi matrix, where 7 = {7?}j = i, k = {kj}1j =1 ; 
7j = defines the local stability of the fixed point along j-th axis, kj = ^s\j 
gives the corresponding oscillations frequency that relates to the wave vector 
magnitude. 



A —A 



D 



d 2 



ef 



dx 2 



X=( Xn.X — .Xj-\ 



g f 1 ^ 9D ef 

dx y 2D 2 ef dx 



From the linear stability analysis it follows that relations for imaginary parts 
k\ = —k2, ks = —k^ and real parts 71 = —72, 73 = —74 are hold. Solv- 
ing Eq.(21) and using these symmetrical properties, we present only positive 
values of k and 7 as functions of the noise intensity a 2 in Fig. 3. Obviously, 




Fig. 3. Positive values of imaginary k (a) and real 7 (b) parts of eigenvalues A 
versus the noise intensity cr 2 , other parameters are: e = 0.2, \x = 0.5, a = 0.2, 
k = 1. Solid dashed and dotted lines correspond to stability of the fixed points 
0(0,0,0,0), (#+, 0,0,0) (x-, 0,0,0), respectively. Plots (a) correspond to lateral 
interaction model, plots (b) are related to x 4 -model. 

if one of the real part related to the fixed point of a homogeneous state is 
positive, then the such point is unstable in the r-space. In the case of zero 
values of all real parts for the fixed point with nonzero imaginary ones, the 
fixed point is a center of a manifold type of limit cycle (with a single non- 
trivial kj) or torus (with more than one nonzero imaginary parts). To relate 
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kj to the profile x(r) or the corresponding derivatives we have solved Eq.(20) 
in the one-dimensional case at the appropriate values for the system param- 
eters. The corresponding description will be given below. As Fig. 3 shows at 
small noise intensities with a 2 < all fixed points are unstable, no pat- 
terns can be formed. If the threshold a\_ is crossed, then the stability of the 
point (x-, 0,0,0) is changed. It becomes a center of a set of tori with two 
frequencies k\ and ki displayed as dotted lines in Fig. 3. When the threshold 
a 2 is crossed, then the stability of the solution 0, 0, 0) is changed again, it 
becomes an unstable one. Therefore, in the interval a\_ < a 2 < a 2 the phase 
space is characteized by a single set of tori formed at the center of the point 
(x-, 0,0,0). At a 2 > a 2 the null state, characterized by the point (0,0,0,0) 
acquires a neutral stability, which is realized till a 2 = a^ Q (see solid lines in 
Fig. 3), where 



K 2 



°™ = < + j a - (22) 

Moreover, at a 2 = a\ + the fixed point (x + , 0,0,0) becomes a center of an 
additional set of tori (dashed lines). Hence, in the interval of the noise intensity 
o"t + < a 2 < a 2 two set of tori are observed. At a 2 = a 2 the stability of 
the point (x + , 0,0,0) is changed, it becomes unstable due to the one of the 
corresponding real part is positive. Therefore, two sets of tori around the 
points O(0, 0, 0, 0) and (x + , 0, 0, 0) are realized till the threshold a 2 . Due to the 
bifurcation related to the noise induced transition oscillations in the vicinity 
of the point (x + , 0,0,0) disappears at a 2 = a 2 and as a result a single set of 
tori around the point (0, 0, 0, 0) can be observed. 

The value a^o as other thresholds except a 2 and a 2 are dependent on the 
parameter k. Above, we have considered the case when a 2 < a^ . From our 
analysis it follows that at small noise intensities or without noise the homoge- 
neous states are absolutely unstable. With an increase in the noise intensity 
two-period solutions of Eq.(20) are realized. At large noise intensity the sys- 
tem is characterized by the unstable homogeneous solution. Therefore, we get 
a reentrant picture of self-organization with the noise growth. The correspond- 
ing phase diagram illustrating patterns formation is shown in Fig. 4. Here we 
plot the parameter a versus noise intensity a 2 at fixed values for other control 
parameters. Unstable states are outside the filled domain. In the filled part 
we denote three sub-domains (a), (b) and (c). In the sub-domains (a) and 
(c) only one set of tori is formed around points (x_, 0,0,0) and (0,0,0,0), 
respectively. The domain (b) is characterized by two sets of tori in the vicinity 
of both points (0,0,0,0) and (x + , 0,0,0). All calculations are performed for 
two choices lateral interactions and x 4 -model of phase separation. Comparing 
Fig. 3a and Fig. 3b, one can see that there are no qualitative changes in the 
dependencies of the eigenvalues of the Jacobi matrix versus the noise inten- 
sity. The situation is similar for these two models. The corresponding phase 
diagrams for two models of phase separations are topologically identical. In 
the x 4 -model the size of the domain of patterns formation in the vicinity of 
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a) & b) ° 

Fig. 4. Phase diagram of reentrant behaviour of self-organization process at e — 0.2, 
\i = —0.5, (3 = 1, K = 1.0. Domains of tori formation are filed: (a) corresponds to a 
tori set around the point (x_, 0, 0, 0); (c) is addressed to a tori set around the fixed 
point (0,0,0,0); double filed domain (b) relates to two set of tori around points 
(x_, 0, 0, 0) and (x + , 0, 0, 0). Plot (a) correspond to lateral interaction model, (b) is 
related to x 4 -model. 

the point 0,0,0) is reduced. 

Let us relate every root hi to the corresponding solution of the stationary 
equation (20), namely to x(r), y = Vx(r), z = V 2 x(r), u = V 3 x(r). To make 
a more accurate analysis we use conditions that a ratio of two different wave 
vectors ki/kj gives a rational number and appropriate real values equal zero, 
ji = 7j = 0, where i ^ j. We plot all solutions in Fig. 5 in a vicinity of every 
fixed point (x_ ,0,0,0), (x + ,0,0,0) and (0, 0, 0, 0) [M It is seen that the profile 
x{r) has one large period related to small values of k. Two-period solutions are 
related to profiles of spatial derivatives of the field x. Therefore, the stationary 
noise sustained structures x st (r) has one period defined by branches with small 
k values in Fig. 3. The corresponding solutions for x 4 -model are of the same 
kind, no qualitatively changes are observed. 

Solutions of the problem (20) with irrational ratio ki/kj are shown in Fig. 6. 
Here we choose the noise intensity values to be about the well known ir- 
rational number with a good precision, ki/kj ~ tt. It follows that in the 
three-dimensional functional space solutions of the problem (20) is a three- 
dimensional line which lies on the torus. This line is an unclosed and the 
corresponding torus is dense filled (the corresponding fractal correlation di- 

1 To identify that the ratio ki/kj is about rational or irrational number we use 
an additional criteria. Due to in numerical calculations all numbers are understood 
as rational ones with given precision in our analysis we have used an algorithm 
allowing to determine a fractal correlation dimension D c of a manifold (torus) in 
the space (x, Vx, V 2 x, V 3 x). In our calculations the value D c ~ 1 corresponds 
to the rational ratio ki/kj (the solution of the problem (20) is a three-dimensional 
closed line (torus) in the space (x, Vx, V 2 x, V 3 x)), whereas D c ~ 2 relates to the 
irrational one (the corresponding solution is an unclosed line on the torus in the 
four-dimensional space) . 
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a) b) c) 

Fig. 5. Solutions of the stationary equation (20) spatial profiles for x(r), y(r) = Vx, 
z(r) = Vy, and u(r) = Vz: (a) a 2 = 0.75832; , x(0) = -0.083, k-/kj = 4; b) 
a 2 = 1.15896, x(0) = 0.425, kjkj = 2; c) a 2 = 1.45(45), x(0) = 10" 5 , ki/kj = 3. 
Other parameters are: a = 0.2, e = 0.2, \x = —0.5, initial conditions for derivatives 
are y(0) = z(0) = u(0) = 0. 




Fig. 6. Solutions of the stationary equation (20) (set of tori) in a space (x, y = Vx, 
z = Vy) at ratio ki/kj ~ tt: a) around fixed point x_ at a 2 = 0.635213; b) around 
fixed point x+ at a 2 = 1.2532123; c) around fixed point xq at a 2 = 1.41767865. 
Other parameters are: a = 0.2, e = 0.2, \x = —0.5, initial conditions for derivatives 
are: y(0) = z(0) = u(0) = 0. 



mension is D c ~ 2.0 ± 0.02). 
4.3 Simulations 

To investigate spatial structures numerically let us consider solutions of the 
Langevin equation (5) numerically. In order to perform numerical analysis, 
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we redefine the model considering a regular two-dimensional lattice with N 2 
points and a lattice spacing i. Then, the partial differential equation (5) is 
reduced to a set of usual differential equations written for an every cell i on a 
grid in the form 

dr BF 1 

= /(*«) + (V^Z^XV^— + -j^=at), (23) 



D(xi) 

where index % labels cells, % = 1, . . . ,N 2 ; the discrete left and right operators 
are introduced 

1 1 1 (24) 

(Vi,)ij = -(Vfl)ji, (VL)y (Vfl)ji = A u = p(6ij+i - 25ij + 6ij-i), 

all sums are taken according to the Einstein rule. A derivative of the free 
energy F is of the form 

For stochastic sources the discrete correlator is of the form (£i(t)£j(t)) = 
2£ 2 a 2 5ij5(t — t'). In the following analysis we use the Stratonovich interpre- 
tation of the Langevin equations (5). In our integration procedure we choose 
i = 1. After discretization scheme described above, the resulting set of cou- 
pled stochastic ordinary differential equations has been integrated numerically. 
To simulate evolution of the system we have used a first-order Euler algo- 
rithm with a time step St = 0.005. Initial conditions for the field were chosen 
from a random Gaussian distribution with fixed average (x(r, 0)) (in vicin- 
ity of the fixed point defined from analysis of the effective potential U e f(x), 
see Fig. 2) and small dispersion, ((Sx(r), 0) 2 ) = 0.003. Typical patterns ap- 
peared in a course of evolution at different values of the noise intensity a 2 are 
shown in Fig. 7. Evolution of the two first moments is shown in Fig. 8. From 
numerical solutions it follows that at small noise intensity (see Fig. 7a) the 
system evolves by a scenario when inside a matrix phase a new phase appears. 
The stationary value of the first statistical moment (x(r,t)) goes from initial 
one (x(r, 0)) ~ x_ to a positive fixed value (see Fig. 8a ). It means that at 
small noise intensity the effective potential U e f(x) has a broken symmetry: 
U e f(x-) > U e f(x + ). It leads to the effect when in a course of evolution the 
value (x) goes to a positive constant. At a 2 = a 2 one has U e f(0) = U e f(x + ). 
Then, starting from the configuration (x(r, 0)) = x_ > with small dis- 
persion a picture kind of phase separation should be realized. Indeed, de- 
spite the picture shown in Fig. 7b is noisy with domains of two phases and 
diffuse interfaces, the statistical average (x) does not change its value, i.e. 
(x(r, 0)) = (x(r,t — > oo)) = x_ (see Fig. 8b). Therefore, here we have a situa- 
tion when dynamics is conserved: / x(r, t)dr = const, but there are no domain 
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Fig. 7. Spatial patterns appeared in a course of evolution at t = 4000: a) a 2 = 0.2; 
b) a 2 = a\ = 1.2757; c) a 2 = 2.0. Other parameters are: e = 0.2, fx = 0.5, k = 1.0, 
P = 1.0, a = 0.2 




Fig. 8. Evolution of first and second statistical moments at a 2 = 0.2 (a), 
g 2 = <Tq = 1.2757 (b) and a 2 = 2.0 (c). Other parameters are: e = 0.2, \i = 0.5, 
K = 1.0, /? = 1.0, a = 0.2 




Fig. 9. An averaged structure function at different values of the noise intensity at 
t = 4000 (the model of lateral interactions): triangles correspond to a 2 — 0.25; 
circles correspond to a 2 = o"q; stars relate to a 2 = 1.5. Other parameters are: 
e = 0.2, /u = 0.5, k = 1.0, (3 = 1.0, a = 0.2 

size growth law as it realizes in phase separation processes. At a 2 > a 2 c (see 
Fig. 7c) when the effective potential has only minimum placed at x = the 
numerical solution (see Fig. 7c) shows massive structures with (x(r,t)) ~ 0.0. 
Despite the system has trivial value for the first statistical moment, the second 
moment is not trivial (see Fig. 8). We calculate the structure function at differ- 
ent values a 2 . The corresponding dependencies are shown in Fig. 9. It is seen 
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that S(k) has one peak located at values k determined by the linear stability 
analysis (see Fig. 3a) and depends on the noise intensity. With an increase in 
a 2 the peak smears, it means that patterns have more diffuse interface. 

The bifurcation values for a 2 illustrating a qualitative change of the station- 
ary distribution (see Fig. 2a) can be calculated numerically. Due to we can not 
investigate a form of the total stationary distribution functional "P^fx] numer- 
ically, our analysis were performed for the quantity 5 = (x) — x mp , where x mp 
is the most probable value related to the position of the maximum for the 
distribution function P st . In such a case 5 plays a role of a criterium for the 
quantitative change of the stationary distribution form. Indeed, considering 
the interval a 2 G [0,af), one choose x mp = that yields 5 = (x) > 0. It 
means that the stationary distribution has a broken symmetry with a large 
contribution related to values x > 0. At a 2 = a 2 the stationary distribution 
has a double degenerated point located at X — X q — X — • It leads to the fact 
that a dependence 5{a 2 ) breaks at a 2 = a 2 and decreases with a 2 growth; 
one has to take into account the value x mp = x_ > in the parameter 5 
calculation at a 2 G [a 2 , a 2 ). At a 2 = <7q one has (x) = x_ and 5 = 0. It 
means that the stationary distribution is of a symmetrical form with respect 
to the central extremum x = x_. With further increase in the noise intensity 
the parameter 5 becomes negative and changes its value abruptly at a 2 = a 2 
where x mp = 0. The corresponding dependencies (analytical and numerical) 
are shown in Fig. 10. It is seen a good correspondence of analytical results with 
computer simulations. 

Previously, we have studied stochastic dynamics when the system evolves into 
the stationary state, described by the functional V s t[x] oc exp(— U e f[x]/a 2 ). 
Let us consider the case when the system relaxes according to the stationary 
distribution. To generate the field x according to this distribution we can use 
the algorithm known variously as the "Langevin method" [31]. The principle 
idea lies in a supposition that the distribution is realized, already. Then, the 
stochastic field x can be obtained as solution of an effective Langevin equa- 
tion d t x = —5U e f[x]/5x + £(r, £), where £(r, £) is a white additive noise with 
intensity a 2 . Solutions of the deterministic equation d t x = —5U e f[x]/5x at 
t — > oo allows to give the stationary patterns at fixed noise intensities. Here 
one assumes that the system is distributed according to an effective probabil- 
ity density functional V e f — Hi (Si(8U e f[x]/8x)), where average is taken over 
initial conditions, % relates to minima positions of the function U e f(x). The 
corresponding patterns at different a 2 are shown in Fig. 11. It is seen that 
at small noise intensity (Fig. 11a) nuclei are formed, at intermediate values 
(a = ctq) a picture type of spinodal decomposition is realized (see Fig. lib), at 
large a 2 (see Fig. 11c) linear defects of dislocations kind are realized. Obtained 
pictures correspond to profiles of the concentration x st , shown in Fig. 5: at 
small noise intensities the pointed profile is observed, whereas at intermediate 
and large a 2 a symmetrical form of profiles are realized. 
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2 
(7 

Fig. 10. The effective order parameter 8 versus noise intensity: solid line corresponds 
to analytical investigation; triangles correspond to simulations on the lattice with 
linear size N = 96; circles are related to N = 192. Other parameters are: a = 0.2, 
e = 0.2, fi = -0.5, f3 = 1.0, « = 1. 




Fig. 11. Stationary patterns at a 2 = 0.2 (a), a 2 = Gq (b) and a 2 = 2(c). Other 
parameters are: a = 0.2, e = 0.2, \x = —0.5, (3 = 1.0, k = 1. 



5 Conclusions 

We have studied a simple deterministic model of reaction-diffusion systems 
which can qualitatively describe stationary noise patterns. In order to study 
noise induced mechanism for pattern formation we add a multiplicative noise 
in a way to satisfy a fluctuation-dissipation relation. We have found that 
the stationary distribution can be obtained exactly. Comparing noise induced 
transitions picture and pattern formation processes it was shown that the 
system follows the entropy driven mechanism by analogy with entropy driven 
phase transitions theory. 

We have explored early stages of the system evolution and found that at 
small times the noise leads to instability of the null state. Considering the sta- 
tionary case we solve the variation problem when the stationary distribution 
functional V st [x] is maximized. Solutions of the problem show that station- 
ary patterns are formed in the vicinity of the maxima of the corresponding 
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stationary density function, P st (x). In linear stability analysis in r-space in 
the stationary case we have found that at small noise intensities the system is 
characterized by unstable homogeneous solutions related to the maxima posi- 
tions of the function P st (x). With an increase in the noise strength the spatial 
structures x st (r) are formed in the vicinity of the local maxima of P s t(x). Ana- 
lytically and numerically we have obtained that at large noise intensity after a 
bifurcation point (the stationary probability density function P s t(x) is of one 
maximum) the noise can sustain stationary patterns. The process of the pat- 
tern formation is well described by the formalism of noise induced transitions 
in zero-dimensional systems, whereas the obtained effect is a some generaliza- 
tion of the noise induced transitions for extended systems. Our study shows 
that at small noise intensity the system manifest a nucleation regime, at fixed 
value of the noise strength a spinodal decomposition is realized, at large noise 
the system exhibits strip patterns with liner defects. Strip structures exist 
in the fixed interval of the noise intensity — large fluctuations destroy the 
patterns. 

Our results are in a good correspondence with experimental data obtained for 
patterns formation in polymer gels with chemical reactions (see Ref. [32]). It 
was shown that in such the reaction-diffusion system patterns can be formed 
only inside the fixed temperature interval (10°C < T < 60°C). Author dis- 
cusses diffusion-induced pattern formation mechanism of volume phase tran- 
sition with pattern formation. He shows that straight strip patterns tend 
to appear in high temperature regime (T « 57°C), and random patterns 
appear at low temperature (T « 30°C). Moreover, the size of patterns de- 
crease in a power law form with an increase in the temperature, L oc T - ^/ 2 , 
j3 = 1.01 ±0.03. In our calculations we have not used the corresponding chem- 
ical reactions, the local dynamics were governed by the Schlogl model [30], we 
have used the model of non-Fickian transport with the field dependent diffu- 
sion coefficient. The noise source satisfies the fluctuation-dissipation theorem 
and has intensity related to the bath temperature. To compare our results 
with above experimental investigations we need to stress that an actual noise 
intensity values are a 2 £ (cr 2 ,<7y ]. We have shown that random stationary 
patterns are realized if a 2 < a 2 < a 2 (experimental photograph in Ref. [32] is 
the same as we show in Fig. lib). At a 2 < a^ the strip patterns are realized 
(see Fig. 11c), the same picture was obtained in Ref. [32]. Moreover, we have 
calculated the dependence L oc (cr 2 ) - ^/ 2 , and obtained j3 = 3.2 ± 0.02. The 
qualitative correspondence in the wavelength dependence is observed. The dif- 
ference in the exponent j3 values can be explained by the fact that we have 
used another kind of the local dynamics model, non-Fickian transport and 
white noise assumptions. 

Obtained results can be applied to study patterns in adsorption/desorption 
processes in metals [29] deposition of a monolayer of molecules [33] and in 
processes of microstructure transformations of materials subject to intensive 
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irradiation. One can await that formation of point defect clusters and clus- 
ters of particles in the system kind of "vacanies+particles" where a phase 
decompositions and "chemical transformations" caused by irradiation can be 
observed. 
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